This book presents the latest numerical solutions to initial value problems and boundary value problems described by ODEs and PDEs. The author offers practical methods that can be adapted to solve wide ranges of problems and illustrates them in the increasingly popular open source computer language R, allowing integration with more statistically based methods. The book begins with standard techniques, followed by an overview of 'high resolution' flux limiters and WENO to solve problems with solutions exhibiting high gradient phenomena. Meshless methods using radial basis functions are then discussed in the context of scattered data interpolation and the solution of PDEs on irregular grids. Three detailed case studies demonstrate how numerical methods can be used to tackle very different complex problems. With its focus on practical solutions to real-world problems, this book will be useful to students and practitioners in all areas of science and engineering, especially those using R.

Mathematical modelling of physical and chemical systems is used extensively throughout science, engineering, and applied mathematics. To use mathematical models, one needs solutions to the model equations; this generally requires numerical methods. This book presents numerical methods and associated computer code in Matlab for the solution of a spectrum of models expressed as partial differential equations (PDEs). The authors focus on the method of lines (MOL), a well-established procedure for all major classes of PDEs, where the boundary value partial derivatives are approximated algebraically by finite differences. This reduces the PDEs to ordinary differential equations (ODEs) and makes the computer code easy to understand, implement, and modify. Also, the ODEs (via MOL) can be combined with any other ODEs that are part of the model (so that MOL naturally accommodates ODE/PDE models). This book uniquely includes a detailed line-by-line discussion of computer code related to the associated PDE model.